The Spectral Energy Distribution of Galaxies 
Proceedings lAU Symposium No. 284, 2011 
R.J. Tuffs & C.C.Popescu, eds. 



©2011 International Astronomical Union 
DOl: 00.0000/XOOOOOOOOOOOOOOOX 



Variation in the dust spectral index across 

M33 

Fatemeh S. Tabatabaei^, Jonathan Braine^, Carsten Kramer^, 
Manolis Xilouris^, Mederic Boquien^, Simon Verley^, 

Eva Schinnerer"^, Daniela Calzetti^, Prancoise Combes^, Prank 
Israel^*^, Christian Henkel^^ and The HerM33es Team 

^Max-Planck-Institut fiir Astronomie, Konigsthul 17, 69117-Heidelberg, Germany 

email: taba@mpia.de 

■^Laboratoire d'Astrophysique de Bordeaux, Uiiiversite de Bordeaux, 33271 Floirac, France 

email:braine@obs . u-bordeauxl . f r 

''institute Radioastronomia Milimctrica, 18012 Granada, Spain 
email:kramer@iram . es 

''institute of Astronomy and Astrophysics, National Observatory of Athens, GR-15236 Athens, 

Greece, emaihxilourisOastro .noa.gr 

^ Laboratoire d'Astrophysique de Marseille, UMR 6110 CNRS, Marseille, France 
emaihmederic . boquienOoamp . f r 

^Department Fi'sica Teorica y del Cosmos, Universidad de Granada, Spain 

emaihsimonSugr . es 

'^Max-Planck-Institut fiir Astronomie, Konigsthul 17, 69117-Heidelberg, GermanyGermany 

email: schinner@mpia.de 

^Department of Astronomy, University of Massachusetts, Amherst, MA 01003, USA 
email:calzetti@astro .umass . edu 

^ Observatoire de Paris, 75014 Paris, France 
email: francoise.combes@obspm.fr 

^'^Sterrewacht Leiden, Leiden University, PO Box 9513, 2300 RA, Leiden, The Netherlands 

email: israel@strw.leidenuniv.nl 

^^Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121-Bonn, Germany 

email: chenkel@mpifr.de 

Abstract. Using the Herschel PACS and SPIRE FIR/submm data, we investigate variations in 
the dust spectral index /3 in the nearby spiral galaxy M33 at a linear resolution of 160 pc. We 
use an iteration method in two different approaches, single and two-component modified black 
body models. In both approaches, /3 is higher in the central disk than in the outer disk similar 
to the dust temperature. There is a positive correlation between /3 and Ha as well as with the 
molecular gas traced by CO (2-1). A Monte-Carlo simulation shows that the physical parameters 
are better constrained when using the two-component model. 
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1. Introduction 

Inferring the nature of dust and dust heating requires knowledge about the dust emis- 
sivity, and its power law index. There is some evidence supporting variation of dust 
emissivity with environmental conditions (e.g. Lisenfeld et al. 2000, Paradis et al. 2009). 
The dust spectral index could change depending on grain properties, e.g. geometry, size 
distribution, and chemical composition. Dust grains might be affected by different phys- 
ical processes like shattering, sputtering, grain - grain collision, condensation and co- 
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Figure 1. Distribution of dust spectral index (right) and column density (in /igcm ^, right) 
across M33, based on the single-component approach. 

agulation. Hence, depending on recent star formation activities, dust/gas content and 
metallicity changes across galaxies, variations in the dust spectral index in inner/outer 
disks or arm/inter-arm regions may occur. Herschel observations of the nearest Scd galaxy 
M33, (HerM33es, Kramer et al. 2010, Boquien et al. 2010) ideally enable us to derive 
distribution of /? across M33. 

We use the Herschel PACS & SPIRE data at 160, 250, 350, and 500/im together with 
the Spitzer MIPS data at 70/im to investigate the physical properties of a standard dust 
model (Big grains+vsg+PAH, Kriigel 2003 and references therein). These wavelengths 
spread over a range where the dust emission mainly emerges from the big grains, which 
are believed to be in thermal equilibrium with the interstellar radiation field. 

2. Analysis and Results 

We use the Newton-Raphson iteration method to derive T , /3 and column density, pixel 
by pixel, in two different approaches: 

1) We assume that dust emitting at the SPIRE bands is mainly heated by the interstellar 
radiation field (ISRF) and obeys single - component gray body radiation. In this case we 
use the data at longer wavelengths (250, 350 and 500/im). 

2) We consider another dust component heated by young massive stars with higher tem- 
perature (warm dust) besides the dust heated by the ISRF (cold dust), hence, dust emits 
as a two-component gray body. Here, we use the data from 70/im to 500//m wavelengths. 
Using a Monte-Carlo simulation, we further investigate the reliability of each approach 
and estimate the errors and uncertainties. We derive the following results: 

Single Component Approach As shown in Fig. 1, in the central disk (i? < 4 kpc), 
the dust spectral index is higher in the central disk (/3 >1.4) than in the outer disk. 
Similarly, both T and column density are higher in the central disk. The flocculent spiral 
structure of M33 is evident in the column density map. The average T and /3 for the 
entire galaxy are 22±9 K and 1.3±0.5, respectively. 




Figure 2. Correlations between tlie dust spectral index and Ha, C0(2-l), and HI integrated 
intensities. The WLS weighted least square fit (red line) and its 5-sigma confidence levels (green 
lines) are shown. The blue lines show the ordinary least square fit (OLS). 

Two Component Approach The cold and warm dust components are separated as- 
suming a variable j3 for the dominant colder component. Although we have to use a 
fixed value of /3warm for the warm dust, we repeat the analysis for different possibili- 
ties of /3warm=l, 1-5, and 2. Interestingly, we find that /3 is again higher in the central 
disk and even shows the arm/inter-arm contrast, irrespective of /3warm assumed. Also 
in all three cases, the cold dust follows spiral arms and dust lanes and is denser where 
the temperature is lower. The warm dust (both temperature and column density) fol- 
lows the distribution of Ha and star forming regions. For the entire galaxy, we derive 
Tcoid= 16 ±3 K, Twarm=40± 13 K, and (5= 1.7 ±0.4 assuming /3^^,ai.,„=2. Our simulations 
show that the two-component approach, generally, constrains the parameters better than 
the single-component approach, particularly the temperature. 



3. Correlation with ISM Tracers 

The M33's outer/inner disk and even the spiral arms are visible in the /3-maps. Having 
similar distribution as of the molecular gas and radiation field traced by Ha, /3 seems to 
have a meaningful change across the galaxy. Figure 2 shows correlations between /3 and 
Ha, C0(2-l) and HI integrated intensities (Gratier et al. 2010) as tracers of the ionized, 
molecular, and atomic phases of the ISM, respectively. In both approaches, /? is found 
to be better correlated with both Ha and CO than with HI. This could indicate that /? 
is not the same in star forming and non-starforming regions in a galaxy. Star formation 
could in principle influence the grain composition (via a change in metallicity) and size 
distribution (via shattering of dust grains due to strong shocks), which could modify the 
dust emissivity index i.e. /3. 
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